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ABSTRACT 

Most  approaches  to  data  association  in  target  tracking  use  a  likelihood-ratio  based  score  for  measurement-to-track  and 
track-to-track  matching.  The  classical  approach  uses  a  likelihood  ratio  based  on  kinematic  data.  Feature-aided  tracking 
uses  non-kinematic  data  to  produce  an  “auxiliary  score”  that  augments  the  kinematic  score.  This  paper  develops  a  non- 
kinematic  likelihood  ratio  score  based  on  statistical  models  for  the  signal-to-noise  (SNR)  and  radar  cross  section  (RCS) 
for  use  in  narrowband  radar  tracking.  The  formulation  requires  an  estimate  of  the  target  mean  RCS,  and  a  key  challenge  is 
the  tracking  of  the  mean  RCS  through  significant  “jumps”  due  to  aspect  dependencies.  A  novel  multiple  model  approach 
is  used  track  through  the  RCS  jumps.  Three  solution  are  developed:  one  based  on  an  a -filter,  a  second  based  on  the 
median  filter,  and  the  third  based  on  an  IMM  filter  with  a  median  pre-filter.  Simulation  results  are  presented  that  show  the 
effectiveness  of  the  multiple  model  approach  for  tracking  through  RCS  transitions  due  to  aspect-angle  changes. 
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1.  INTRODUCTION 

Feature-aided  tracking  (FAT)  is  a  technique  where  feature  data  are  used  to  enhance  the  data  association  process  in  tracking. 
In  the  context  of  track-to-track  correlation,  the  use  of  features  is  called  feature-aided  track  correlation  (FATC).  Kinematic 
scores,  based  on  distance  metrics,  are  computed  in  the  classical  approach  for  determining  the  likelihood  that  a  particular 
report  should  be  associated  to  a  given  track.  When  feature  data  (e.g.,  RCS,  intensity,  length,  etc.)  are  provided  with  the 
report,  it  is  possible  to  compare  it  to  a  current  estimate  of  that  feature  for  the  given  track.  Thus,  association  is  enhanced 
when  the  report  feature  is  similar  to  the  track  feature,  and  discouraged  when  it  is  different.  The  main  challenge  in  FAT 
or  FATC  is  to  formulate  a  “score”  (typically  a  likelihood  ratio)  that  can  be  used  in  the  association  algorithm  to  quantify 
the  similarity  during  association  processing.  Further,  as  the  feature  associated  with  a  track  changes  over  time,  one  needs  a 
mechanism  to  “track”  the  changes.  In  the  case  of  the  RCS  feature,  the  implementation  of  a  dynamical  model  is  difficult  or 
prohibitive  because  of  a  lack  of  information  about  the  target  type  and  aspect-angle  dependency. 

A  number  of  different  approaches  to  FAT  have  been  developed  in  the  literature.  Drummond1-2  describes  a  likelihood 
function  for  use  in  data  association  that  incorporates  the  feature;  errors  in  measuring  the  feature  are  assumed  to  be  Gaussian 
(thus  the  likelihood  involves  a  Chi-square  parameter),  and  the  feature  mean  is  assumed  to  constant  over  time.  Lerro  and 
Bar-Shalom3  integrated  an  amplitude  likelihood  ratio  into  the  standard  PDAF  association  probabilities,  which  affected  the 
association  probabilities  by  favoring  a  high-amplitude  measurement.  Singer  and  Coursey4  developed  a  non-parametric 
FAT  approach;  this  is  a  useful  approach  when  nothing  is  known  about  the  feature  measurement  density  function.  Agate 
and  Sullivan5  provide  a  general  derivation  of  measurement-to-track  association  hypotheses  that  include  features  and  kine¬ 
matic  data;  the  work  assumes  the  availability  of  a  library  of  class-dependent  feature  measurement  PDFs  based  on  target 
classes  and  aspect  angles  (the  main  application  envisioned  is  high  range  resolution  radar).  Nguyen*’  et  al.  developed  an 
approach  that  accumulates  target  classification  information  obtained  from  high  range  resolution  radar  to  improve  associa¬ 
tion  scores.  Lancaster,  Blackman,  and  Taniguchi '  developed  a  technique  using  Dempster-Shafer  methods  for  incorporating 
target  classification  data  into  the  data  association  score. 

In  this  paper,  we  address  the  specific  problem  of  using  radar  signal  features  (SNR  and  RCS  measurements  from  a 
narrowband  radar)  for  augmenting  the  track  score  used  in  the  data  association  problem.  There  are  two  main  components 
to  the  work:  (i)  the  signal  score  likelihood  ratio  formulation,  and  (ii)  the  RCS  tracking  technique.  The  general  signal 
score  likelihood  ratio  for  data  association  used  to  augment  the  kinematic  score  is  described  in  Blackman  and  Popoli.8 
Likelihood  ratios  for  specific  RCS  fluctuation  models  are  derived  following  the  work  of  Blackman  and  Fitzgerald.9  For 
the  RCS  tracking  technique,  we  advance  the  multi-model  RCS  estimation  approach  in  Ref.  9  that  used  a  smoothing  filter. 
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We  extend  the  work  in  two  ways.  First,  we  examine  a  variation  on  the  approach  by  using  median  filters  instead  of  the 
smoothing  filter.  Second,  we  formulate  a  new  switching  model  approach  that  uses  an  IMM  filter  combined  with  a  median 
filter.  The  median  filter  is  ideal  for  processing  data  that  has  a  “large  tail”  density  (e.g.,  Swerling  I). 

This  paper  is  organized  as  follows.  In  Section  2,  we  introduce  some  notation  used  throughout  the  paper.  Section  3 
develops  the  signal  score  for  three  different  RCS  fluctuation  models:  Swerling  I,  Swerling  III,  and  log-normal.  Section  4 
develops  three  multi-model  approaches  for  RCS  estimation.  Section  5  presents  simulation  results  that  compare  the  RCS 
estimation  approaches.  Finally,  Section  6  summarizes  the  work. 

2.  NOTATION 

For  the  convenience  to  the  reader,  we  summarize  here  the  notation  used  in  the  subsequent  developments. 

L^(k)  =  the  log-likelihood  function  for  measurement-*,  track-j,  at  time  tk- 
fj'fr  ( k)  =  the  kinematic  portion  of  the  log-likelihood  function. 

Lg  ( k )  =  the  signal  portion  of  the  log-likelihood  function. 

Rjc  =  the  observed  SNR  of  the  yth  measurement  (linear  units)  at  time  tk- 
XJk  =  the  RCS  of  the  yth  measurement  (in  m2)  at  time  tk- 
crJk  =  the  RCS  of  the  jth  measurement  (in  dBsm)  at  time  tk- 

&k\k_i  =  the  estimate  of  the  target  RCS  associated  with  the  vth  track  (in  dBsm)  at  time  tk  using  data  through 
time  tk- 1. 


3.  SIGNAL  SCORE  MODEL 

The  track  score  (log-likelihood)  for  the  ?'th  track  updated  with  the  jth  measurement  at  time  is  described8, 9  to  be 

Lij  (k)  =  L\k  -  1)  +  A Lij  (k)  (1) 


where 


A  Lij{k) 


ln(l  —  Pd),  no  track  update 

A Llj/(k)  +  Aid?  (fc),  track  update 


(2) 


Here,  AL^(fc)  is  the  incremental  kinematic  portion  of  the  score  update  while  A Ll/l  ( k )  is  the  incremental  signal  portion  of 
the  score  update.  The  expression  for  the  incremental  signal  score  is  given  to  be 


A  Ll/i  ( k )  =  In 


fp{Ri\HllRj>Rth)\ 


(3) 


where  R{  is  defined  to  be  the  observed  SNR  of  the  jth  measurement  at  time  k.  Here,  p(R3k\Hi,  R3k  >  Rth)  is  the 
probability  density  function  (PDF)  for  the  SNR  given  Hi,  the  hypothesis  that  a  signal  is  present,  and  given  that  the 
observed  SNR  exceeds  the  SNR  detection  threshold  Rth-  Also,  p{R3k\Ho,  Kj  >  Rth)  is  the  PDF  for  the  SNR  conditioned 
on  Hq,  the  alternate  hypothesis  of  noise  only  (or  for  clutter),  and  on  the  condition  that  the  SNR  exceeds  the  threshold. 
Since  >  Rth)  =  p{R3k\Hi)/Pd ,  and  p/R/\H(),  sftj.  >  Rth)  =  p{R3k\H0)  /  Pfa,  then  (3)  can  be  simplified  to 

A Lg  (fc)  =  In  (j)(RJk\Hi) /p(3?^,|fTo)^  •  However,  the  term  In  {Pd/ Pfa)  is  often  computed  with  AL^(fc)  when  the  score 
is  implemented  in  software.  Furthermore,  nominal  values  for  Pd  and  Pfa  are  frequently  used.  Thus,  for  this  paper  we  will 
focus  on  the  computation  of  the  “augmenting  signal  score”  that  is  defined  by 


( k )  A  In 


{p(Rl\H1,Rj>Rth)\ 

\p(Ri\H0,Ri>Rth)J 


(4) 


3.1.  Score  for  the  Swerling  I  Fluctuation  Model 

Four  fluctuation  models  for  the  target  RCS  have  been  developed  by  Swerling.10  The  Swerling  I  model  is  a  scan-to-scan 
(slow)  fluctuation  model11  that  is  appropriate  for  targets  with  many  independent  scatting  points  of  which  no  single  or  few 
predominate.  The  Swerling  II  model  is  a  pulse-to-pulse  (fast)  fluctuation  model  with  the  same  RCS  fluctuation  density; 
the  model  is  used  in  cases  where  the  radar  has  a  very  low  PRF,  or  for  targets  where  very  small  changes  in  orientation  cause 
large  RCS  changes. 

The  Swerling  I  density  for  the  RCS  is  given  by 


X^°  (5) 

where  X  is  the  mean  RCS  value.  From  this  density  for  the  RCS,  the  density  for  the  SNR  for  a  Swerling  I  model  (for  single 
pulse  integration)  can  be  shown  to  be  exponentially  distributed,8 

=  “p(rfl)  «s> 

where  ft  is  the  mean  SNR.  Here,  H\  is  used  to  denote  the  hypothesis  that  the  return  is  from  a  target  and  not  noise. 

Given  that  the  return  has  exceeded  the  SNR  threshold  of  ft,,,  to  be  detected,  then 

_  it  (  -SftJ  \ 

v{k\hx&,k  >  »t/o  =  c '  m ' exp  ( rrl )  (7) 


p(Xl\HuX)  =  -  exp 


where  the  normalizing  factor  is  given  by  C  =  p(fflk\H\)d$Vk  =  exp  (— 5ft,/,/(  1  +  5ft))  Thus, 


p(K  //, .  ft.  ft/.  >  ft,/,! 


1 

m 


•  exp 


(-K  + 

^  l  +  K  ) 


(8) 


Meanwhile,  the  PDF  for  the  case  of  noise  only  is  given  by  p(5ft-[|fTo)  =  exp(— 5ft*.).  A  clutter  density  (e.g.,  Weibull) 
could  be  used  instead,  if  appropriate  for  the  application.  Given  that  the  noise  exceeded  the  SNR  threshold  of  ft, h  to  be 
detected,  and  since  exp(5ft)d5ft  =  exp(— 5ft,^),  then  the  SNR  density  for  the  noise  given  it  exceeded  the  threshold  is 

given  by  p(5ft{. | f?o ,  5ft*.  >  5ft,/,)  =  exp(-5ft^.  +  5ft//,). 

Using  this  density  function,  the  augmenting  signal  score  update  becomes 


L3t  (ft)  =  In 


/pi  ft).  //;.ft'  >  ft,/,A 


-ln(l  +  r)  + 


(5ft1  -  5ft^)5ft* 
l  +  5ft* 


(9) 


To  compute  this  score,  one  needs  to  know  5ft®,  the  mean  SNR  for  the  ith  target.  This  can  be  computed  as  5ft*  =  C'y  X ' , 
where  X1  is  the  mean  RCS  and  Co  is  constant  that  is  a  known  function  of  radar  parameters  (power,  wavelength,  range, 
etc.).  In  general,  X1  is  unknown,  but  this  can  be  estimated  as  we  will  discuss  in  Section  4.  The  challenge  is  that  the  RCS 
can  suddenly  change  based  on  the  aspect  angle  to  the  target. 

3.2.  Score  for  the  Swerling  III  Fluctuation  Model 

The  Swerling  III  model  for  target  RCS  is  a  scan-to-scan  (slow)  fluctuation  model11  that  is  appropriate  for  targets  that  have 
one  main  scattering  element  that  predominates  together  with  many  smaller  independent  scattering  elements.  This  model  is 
sometimes  assumed  for  radars  with  lower  frequencies  (e.g.,  below  1.0  GHz)  and  for  small  streamlined  aircraft  targets.  The 
Swerling  IV  model  is  a  pulse-to-pulse  (fast)  fluctuation  model  that  uses  the  same  RCS  density. 


The  PDF  for  the  Swerling  III  RCS  is  given  by 


_  AV  / _ 

p(Xl\H1,X)=J^-exPl^-iry  X>0  (10) 

where  X  is  the  mean  RCS.  The  corresponding  pdf8  of  the  observed  SNR  for  a  Swerling  III  target  (for  single  pulse  integra¬ 
tion)  is  given  by 


1  +  2/kJ 


(TTW'exp 


(irfe) 


4(2 +  £+&&’)  /-  25ftA 

(2  +  5ft)2  6XP  ^2"+! J 


(11) 


where  the  mean  SNR  is  5ft.  Given  that  the  return  was  detected,  then  the  SNR  must  be  greater  than  the  threshold*  The 
normalization  component  is  derived  as 


/-(■ft'  //;.  fth/ft,'.  = 


(5ft2  +  45ft  +  25R5Rt/j  +  4) 

(2  +  !ft)2 


exp 


/-25ft»A 

V  2  +  5?  y 


Therefore,  the  density  for  the  SNR  given  it  exceeds  the  threshold  is 


4(2  +  K  +  5ft5ft£)  _  /  2(ft'  'ft//, !  \ 

(2  +  5ft)(5ft2  +  45ft  +  25ft5ftt,l  +  4)  6XP  ^  2  +  5ft  j 


(12) 


(13) 


To  compute  the  denominator  of  the  log-likelihood  ratio,  we  again  use  the  noise  d4ensity  given  that  the  SNR  is  above  the 
detection  threshold,  p(5ft{|7To,  >  5 Rth)  =  exp(— 5ft)j.  +  5 ftth)- 

Using  these  density  functions,  the  augmenting  signal  score  becomes 


L1!  ( k )  =  In 


//-('ft'  //;. 'ft'. 'ft[  >)(„,)' 

V  p(5fti|JTo,5fti>5ftt/() 


=  In 


4(2  +  5ft4  +  5ft*5ftJfc ) 


(2  +  5ft*)  ((5ft*)2  +  45ft*  +  2RiRth  +  4) 


-  9Efh) 
2  + 5ft* 


(14) 


Again,  to  compute  this  score  one  needs  to  know  5ft*.  An  estimate  can  be  derived  if  one  tracks  the  mean  RCS  of  the  target. 
Sudden  changes  in  RCS  based  on  the  target  aspect  angle  present  the  challenge  in  estimating  X1. 


3.3.  Score  for  the  Log-Normal  Fluctuation  Model 

Another  accepted  RCS  fluctuation  model  is  the  log-normal  distribution,11 

p(Xl\H1,Xx)  = 


2<t2 


In 


(15) 


where  X  is  the  median  value  of  X,  and  ?  is  the  standard  deviation  of  \w(X/ X).  Apparently,  this  model  can  approximate 
the  radar  cross  section  of  some  satellites,  cylinders,  plates,  and  arrays.  With  a  change  of  variables,  cr  =  ln(X),  then  we 
would  find  that  the  PDF  for  cr  is 


P(al\Hi’VX2) 


of 


(16) 


where  a  =  ln(X).  Since  we  typically  report  the  RCS  is  units  of  dBsm,  then  we  can  simply  use  (16)  with  crj,  =  10  log(A5^) 
if  we  define  y  to  be  the  standard  deviation  of  10  log  (A* /A )  instead.  This  model  was  adopted  in  the  report  by  Blackman9 
for  forming  a  signal  score  log-likelihood  ratio. 


The  SNR,  5fti,  and  the  radar  cross  section,  XJk  (m2  units)  or  ai  (dBsm  units),  can  be  related  in  terms  of  the  known  radar 
parameters  Co, 


5ft{)  =  C0  •  XI  =  C0  •  exp 


l  101og(e) 


=  /K) 


(17) 


Note  that  the  “signal  detection  threshold”  is  typically  given  as  7 th  =  5 Rth  •  Pn,  where  Pn  is  the  noise  power. 


Now,  given  an  assumed  probability  density  for  the  RCS  in  (16),  we  can  transform  though  the  function  =  fieri)  which 
is  defined  in  (17)  to  obtain  the  probability  density  function  for  the  SNR, 


Px  (K  I  Hu  a,*)  =Pa  (r\K)  I  Hu<t,<;2 


9f  (y) 


dy 


(18) 


y=K 


Applying  the  transformation,  we  obtain  the  following  density  function  for  the  SNR  of  the  y  th  measurement  (given  a  signal 
is  present). 


p(K\Hi’a^)  =  10 1°$0^  P(al\Hi’a’S2) 


(19) 


where  aJk  =  f  1(3?^,)  =  101og(e)  ln(Jf^,/(7o).  Then,  assuming  the  density  for  the  RCS  from  (16),  the  SNR  density 
simplifies  to 


2<r2 


V^TT^i 

Since  the  detection  is  only  obtained  when  the  return  exceeds  the  threshold  -R/ /, ,  we  must  scale  the  density  by  a  constant. 


(20) 


C  = 


p{Wk\H1:a,<;2)  =  Ml  +  erf 


\/2(cr  -  ath) 

2c 


(21) 


J3tth  2 

Zn  \  \  /  / 

where  <jth  =  101og(e)  ln(5Rt^/Co).  Therefore,  the  pdf  for  the  SNR,  given  the  return  exceeded  the  threshold,  is  given  by 

J  2  ' 

k 

2<;2 


PiKIHud,  c,  >  Vtth)  =  exp  (  {ai  aY 


(22) 


The  null  hypothesis  density,  given  that  the  return  exceeds  the  required  signal-to-noise  is  p(^Rk\H0,  >  With)  = 

exp(— JRj,  +  9tth). 

Inserting  the  density  functions  into  the  likelihood  ratio  gives  us  the  augmenting  signal  score. 


(k)  =  Sftjj.  -  SRtft  -  In ~  ln  C*  +  ln(10  loSio(e))  ~  ln(v^7r?*) 


K  -d‘)2 

2(C)2 


(23) 


Implementation  of  this  solution  requires  the  knowledge  of  C",  <r\  and  C-  The  first  two  come  from  radar  system 
parameters!.  We  note  that  for  most  cases  that  \n(Cl)  «  0,  and  therefore  can  effectively  be  ignored.  The  median  RCS,  <f\ 
has  to  be  estimated  since  it  is  a  target-dependent  parameter.  The  variance  (C)2  is  an  RCS  fluctuation  model  parameter;  it 
could  be  assumed  or  estimated. 


3.4.  Comparisons  of  RCS  and  Median  Filtered  RCS  Fluctuation  Distributions 

Several  of  the  RCS  filtering  techniques  developed  in  Section  4  make  use  of  median  filters.  Thus,  it  is  useful  to  relate  the 
median  value  to  the  mean  value  of  the  Swerling  distributions  discussed  above.  To  find  the  relationship,  we  integrate  the 
Swerling  density  as  follows,  y  =  J0T  p{X\H\,  X)  dX.  Then  the  median  is  2fme(jian  =  T  such  that  y  =  0.5.  Solving  this 
expression  for  the  Swerling  I  RCS  distribution,  we  find  that  Xme<jian  =  1-4427  x  X.  For  the  Swerling  III  distribution,  we 
find  that  Xmedian  =  1.1916  x  X. 

To  understand  the  effect  of  median  filtering  on  RCS  data,  we  performed  simulations  of  the  median  filtered  RCS  dis¬ 
tributions  (window  size  w  =  3)  for  the  Swerling  I  and  Swerling  III  models.  Based  on  a  comparison  of  plots,  the  median 
filtered  Swerling  I  RCS  data  closely  matches  the  distribution  for  the  log-normal  distribution  with  a  standard  deviation  of 
3.0  dBsm.  For  Swerling  III,  the  median  filtered  distribution  closely  matches  the  log-normal  distribution,  but  here  the  stan¬ 
dard  deviation  of  the  log-normal  is  2.0  dBsm  instead.  Thus,  one  can  conclude  that  median-filtered  Swerling  I  and  III  RCS 
data  has  a  distribution  that  approximates  a  log-normal  distribution;  this  characteristic  is  advantageous  when  developing 
RCS  filtering  methods  that  use  a  median  pre-filter. 

^If  the  signal  scoring  is  done  remotely  from  the  radar,  then  it  is  possible  that  the  tracker  may  not  have  access  to  these  parameters. 
The  use  of  approximations  may  be  necessary. 


4.  MULTIPLE  MODEL  RCS  ESTIMATION 

To  compute  three  signal  score  expressions  in  (9),  (14),  and  (23)  requires  and  estimate  for  W  or  equivalently  X1 .  Since 
RCS  is  a  target-dependent  variable,  and  SNR  is  a  function  of  range  and  the  instantaneous  radar  parameters,  it  makes  sense 
that  one  should  track  the  RCS  of  the  target.  From  this  estimate  one  can  infer  the  predicted  SNR.  In  this  section,  we 
describe  three  methods  for  predicting  the  target  mean  RCS  using  tracking  methods.  The  preferred  approach  would  be  to 
have  a  catalog  of  all  target  RCS  patterns  and  to  “look  up”  the  mean  RCS  based  on  a  target  type  estimate  and  an  aspect 
angle  estimate.  This  is  effectively  a  target  classification  problem,  which  has  been  treated  as  an  Hidden  Markov  Model 
(HMM)  problem  in  the  literature.12, 13  However,  maintaining  a  catalog  of  all  target  types  is  often  prohibitive,  therefore 
our  approach  is  to  attempt  to  estimate  the  mean  (or  median)  RCS  based  on  recent  measurements  and  using  no  a  priori  data 
about  the  target.  The  difficulty  in  estimating  the  RCS  is  that  the  mean  RCS  can  rapidly  change  as  the  aspect  angle  to  the 
target  changes.  For  example,  large  returns  can  be  obtained  when  directly  aligned  with  the  nose  section  of  an  aircraft,  while 
smaller  returns  can  be  obtained  at  angles  slightly  off  the  noise.  To  handle  these  “jumps”  in  the  RCS  as  the  aspect  angle 
changes,  three  multiple-model  estimators  have  been  developed.  The  three  approaches  are  compared  in  Section  5. 


4.1.  Two-Model  RCS  Estimator  with  a- Filters 


An  approach  described  in  Ref.  9  uses  a  Bayesian  model  to  allow  for  good  “smoothing”  for  stable  RCS  conditions  and 
“jumps”  when  rapid  RCS  changes  occur  due  to  aspect-angle  dependencies.  Two  estimates  (short-  and  long-term)  are 
obtained  using  two  a-filters.  Define  aJk  to  be  the  RCS  measurement  associated  with  the  yth  measurement  (in  dBsm  units). 
Let  irjufe _x  and  d’^k_1  be  the  short-term  and  long-term  estimates  of  the  RCS^,  respectively,  associated  with  the  ?'th  track. 
Log-normal  probability  densities  are  assumed  for  the  two  RCS  model  hypotheses. 
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Here,  hypothesis  hi  corresponds  to  the  short-term  model  and  hypothesis  h2  corresponds  to  the  long-term  model.  Also,  ?2 
and  ^2  are  variances  associated  with  the  short-  and  long-term  RCS  models.  Using  Bayes’  theorem,  the  a  priori  probabilities 
at  time  tk  are 


Pk(hi)  = 
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where  pk{hn)  is  the  computed  probability  for  filter  mode  hn  using  data  up  to  time  tk .  Initial  prior  mode  probabilities  are 
set  to  be  po(hi)  =  po(h2 )  =  0.5.  Note  that  it  is  possible  that  p(a:k\hi)  — >  0  or  p(aJk\hi)  — >  0,  hence  in  the  calculation 
of  (25)  one  must  place  a  floor  value  pmin  ~  0.01  on  the  mode  probabilities;  otherwise,  the  filter  could  “get  stuck”  in  one 
mode. 


The  short  and  long-term  estimates  of  the  RCS  are  computed  using  an  cr-filter  with  all  the  data  up  to  time  tk, 

K]k  =  (!  -  «iK|fc-i  +  «ioi,  vk\k  =  (!  -  «2 )alk\k-i  +  «2 ai  (26) 

where 

Oil  =  max  (l  -  exp(— Tfc_i/n),  a“m)  ,  a2  =  max  (l  -  exp(-Tfc-_i/r2),  a“ln)  (27) 

and  Tk_i  =  tk  —  tk- 1  is  the  sampling  interval  (time  since  last  update  of  the  /th  track).  The  time  constants  t-[  and  t2 
are  user-defined  parameters  that  should  be  selected  based  on  the  aspect-angle  change  rate  of  the  target  (target  range  and 
velocity)  and  the  sensor  update  rate.  Also,  it  can  be  useful  to  place  upper  and  lower  limits  on  (26)  to  keep  them  within  a 
range  expected  for  targets  of  interest  and  to  mitigate  the  impact  of  outliers. 

Finally,  the  estimated  RCS  is  computed  to  be  the  weighted  average  of  the  short-  and  long-term  estimated  RCS  values. 


°k\k  =  *1) k  '  Pk(h i)  +  yk\k  ■  Pk(h 2) 


(28) 


When  this  blended  RCS  estimate  is  used  to  evaluate  the  likelihood  score,  the  propagated  version  is  used  when  the  next 
frame  of  data  is  received.  However,  without  a  model  for  the  RCS  change  with  time,  the  propagated  RCS  estimate  is  simply 
the  last  updated  value,  dlk+1^k  =  olk |fc. 


*If  we  assume  the  log-normal  distribution,  then  a  is  an  estimate  of  the  median  value  of  the  RCS. 


4.2.  Modified  Two-Model  RCS  Estimator  with  Median  Filters 

A  modification  of  the  approach  shown  in  Section  4. 1  was  developed  by  replacing  the  a-fi  Iters  shown  in  (26)  with  median 
filters.  The  motivation  for  doing  so  was  that  median  filters  are  known  in  the  statistical  community  to  be  a  good  choice  for 
processing  data  with  significant  outliers.  RCS  distributions  such  as  Swerling  I  have  a  significantly  large  “tail,”  thus  there 
is  a  significant  possibility  of  receiving  large  RCS  measurements  relative  to  the  mean  RCS. 

The  filter  is  implemented  by  computing  the  median  value  of  the  RCS  over  the  recent  w i  (short)  or  tv  2  (long)  windows 
of  measurements, 


crk \k  =  median 


jk- 


01+1 


k— iui+1  ’  ' 


'k\k 


=  median 


(a 


Jk—LU  2  +  1 
k  — 1^2  +  1  5  '  *  ’  5 


(29) 


These  estimators  replace  the  a-filters  shown  in  (26).  Notice  here  that  we  have  augmented  the  superscript  notation  such 
that  akk  now  specifies  some  particular  measurement  jp-  on  frame  k ,  and  this  measurement  index  may  be  different  from  the 
measurement  index  jk-i  on  frame  k  —  1. 

The  output  estimate  of  the  target  RCS  is  computed  using  the  Bayesian  blended  RCS  as  shown  in  (28).  Further,  the 
predicted  RCS  on  frame  k  +  1  is  also  defined  by  crlk+1\k  =  &k\k  since  no  RCS  dynamics  model  is  available. 


4.3.  IMM  Filter  RCS  Estimator 

We  present  now  a  new  approach  to  the  multiple  model  RCS  estimation  problem  based  on  an  Interacting  Multiple  Model 
(IMM)  filter.  The  IMM  technique  is  well  known  to  be  the  leading  choice  in  target  tracking  where  multiple  “maneuver 
modes”  are  present.  In  the  RCS  estimation  context,  a  “jump”  in  the  RCS  value  constitutes  a  new  mode  of  the  RCS  state. 

4.3.1.  Jump-Bias  Model  IMM  Approach 

The  general  idea  of  the  IMM  application  for  RCS  tracking  is  to  implement  a  set  of  Kalman  filters  with  bias  inputs  that 
correspond  to  possible  jump  values.  Thus,  when  a  large  RCS  step  is  experienced,  the  filter  will  “switch”  to  the  mode  with 
corresponding  bias  that  best  matches  the  step.  Two  kinds  of  filter  models  can  be  used:  nearly-constant  position  (NCP),  and 
nearly-constant  velocity  (NCV).  The  former  is  useful  when  the  mean  RCS  value  not  changing  with  time,  and  the  later  is 
useful  when  the  mean  RCS  is  slowly  varying  (rapid  changes  are  handled  by  the  multiple  models). 

The  ith  model  state  equation  for  the  RCS  estimation  problem  is  defined  as 

4+i  =  Ffcx^  +  b'  +  wfc  (30) 

where  x'k  and  F;,:  are  discussed  below  for  the  cases  of  NCP  and  NCV,  and  w k  ~  A/"(0,  Qfe)  is  the  process  noise.  The  “jump” 
bias  b'  associated  with  the  / 1 h  model  is  either  a  vector  or  scalar  depending  on  whether  it  is  an  NCV  or  NCP  model.  In  the 
case  of  the  NCV,  then  it  is  a  vector  b'  =  \b‘  ,  0]  ,  where  //  is  the  RCS-mean  jump  value  specified  for  the  ith  model.  In  the 
case  of  NCP,  then  the  bias  is  the  mean  RCS  scalar,  b1  =  bl.  By  defining  a  family  of  jumps  {61,  b2, . . . ,  bNj,  with  bn  =  0 
for  some  n  £  {1, . . . ,  N },  then  the  IMM  allows  the  filter  to  rapidly  transition  to  RCS  steps  whenever  1  /  n  is  the  best 
model,  and  while  providing  smoothing  whenever  the  model  i  —  n  is  the  best  model.  Note  that  the  RCS  jumps  need  not 
match  the  b1  values  since  blending  of  states  takes  place  and  “in  between"  values  will  result  from  the  combination  of  two 
states. 


NCV  Model.  The  nearly-constant  velocity  model  has  a  state  that  includes  RCS  and  RCS-rate,  xp-  =  \rrk,  bp]r .  Thus, 
the  estimate  of  the  mean  RCS  is  extracted  as  the  first  component  of  the  state.  The  RCS-rate  can  be  used  for  predicting  the 
RCS  at  a  future  time,  if  necessary.  The  state  transition  matrix  and  the  process  noise  matrix  are  given  by 
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where  Tp.  =  tk+ 1  —  ffc  and  qo  is  a  user-defined  process  noise  intensity  parameter. 


NCP  Model.  The  nearly-constant  position  model  has  a  state  that  includes  only  the  RCS,  Xfc  =  CTk,  thus,  the  state  is  a 
scalar.  The  state  transition,  Ffc  =  1,  is  also  a  scalar.  Likewise,  the  process  noise,  vv/.  =  r/o 7 ,  is  a  scalar.  Hence,  the  RCS 
state  equation  (30)  for  the  NCP  model  can  be  written  as 


Ofc+i  —  Ofc  +  bl  +  Wk  (32) 

Note  that  the  level  of  “smoothing”  versus  “adaptation  to  change”  in  this  model  can  be  controlled  by  the  choice  of  the  value 
of  q0. 


IMM  Implementation.  The  blended  IMM  output  RCS  state  estimate  is  computed  as 

N 

^  (33) 

i=l 

where  fi\  is  / th  mode  probability,  and  where  xk\k  is  the  Kalman  filter  state  estimate  based  on  the  ith  model.  The  RCS 
estimate  &k\ k  is  then  extracted  from  xk\k.  If  the  IMM  configuration  includes  an  NCV  filter,  then  a  dynamics  model  for  the 
RCS  is  employed  and  a  predicted  RCS  value  <xfc+1|fc  could  be  generated  from  X/i:+1  =  F /,:X/;:|fc.  Other  details  on  the  IMM 

filter  algorithm  implementation  are  well  known  and  can  be  found  in  Ref.  14. 

4.3.2.  Median  Pre-Filter  for  the  IMM  Implementation 

A  unique  aspect  of  the  IMM  filter  implementation  for  RCS  estimation  is  the  application  of  a  median  pre-filter.  The  reason 
for  the  pre-filtering  is  that  RCS  fluctuation  noise,  with  the  large  tail  distribution,  can  degrade  the  performance  of  the 
Kalman-based  IMM  filter.  By  applying  the  median  pre-filter,  outliers  are  suppressed  and  the  error  distribution  becomes 
more  log-normal-like.  The  median  pre-filter  was  successfully  used  in  a  glint  noise  filtering  problem.15  We  have  not 
applied  the  velocity  prediction/correction  idea  developed  in  Ref.  15  because  it  is  not  really  needed  in  this  application;  the 
RCS-rate  is  typically  not  large  and  median  filter  window  is  short  (e.g.,  w  =  3). 

The  pre-filter  processing  is  implemented  as  follows.  If  akk  is  the  jfcth  measurement  at  time  ffc,  then  the  output  mea¬ 
surement  is  formed  as 

aJk  =  median  (af  , , . . , (34) 

where  w  is  the  length  of  the  median  filter  window.  The  output  di  is  passed  to  the  IMM  filter  as  the  RCS  measurement  input 
at  time  ffc,  unless  the  measurements  is  exactly  the  same  as  dk_l  (which  can  happen  because  the  same  measurement  can  be 
the  “median  value”  on  consecutive  frames).  When  the  output  measurement  is  identical,  then  no  measurement  is  passed  to 
the  IMM  filter  (i.e.,  we  coast  the  filter).  This  prevents  correlation  due  to  repeated  measurements  from  building  up  in  the 
IMM  filter  estimates. 

It  is  important  to  note  that  the  median  pre-filter  has  the  effect  of  changing  the  input  distribution  such  that  the  filtered 
RCS  estimates  <5fc|fc  converge  to  the  median  RCS  value.  Thus,  if  one  requires  the  mean  value  for  computation  of  the 
likelihood  ratio  or  for  RMS  accuracy  evaluations,  the  median  estimate  must  be  converted  to  a  mean  estimate  as  discussed 
in  Section  3.4. 


5.  SIMULATION  RESULTS 

To  evaluate  the  three  RCS  tracking  filters,  we  constructed  two  simulation  cases  that  test  the  ability  of  the  filters  to  transition 
overjumps  and  accurately  estimate  the  RCS. 

5.1.  Scenario  Descriptions  and  Filter  Parameters 

Two  simulation  cases  were  constructed.  Case  1  has  static  RCS  values  except  for  jumps  of  10  dBsm  and  20  dBsm;  the  case 
will  show  the  ability  of  the  filters  to  transition  quickly  to  new  RCS  values  while  still  providing  smoothing  during  constant 
RCS  periods.  Case  2  consists  of  slowly  varying  RCS  values  with  some  large  jumps.  This  case  will  show  the  ability  of  the 
filters  to  track  the  RCS  and  respond  to  jumps.  RCS  samples  were  generated  using  a  Swerling  I  distribution,  with  an  update 
rate  of  1.0  Hz. 


Table  1.  Filter  parameters  for  simulation  study. 


Filter  Parameter 

Value 

Q-Filter 

Time  constant  n 

3.0  sec 

Time  constant  t2 

10.0  sec 

Model  SDV  cn 

5.0  dBsm 

Model  SDV  a2 

5.0  dBsm 

Min  Mode  Probability  Pmin 

0.01 

Median  Filter 

Short  Window  wi 

3  samples 

Long  Window  w2 

11  samples 

Model  SDV  cn 

3.0  dBsm 

Model  SDV  a2 

3.0  dBsm 

Min  Mode  Probability  Pmin 

0.01 

IMM  Filter 

Median  Pre-Filter  Window  w 

3  samples 

RCS  Measurement  Covariance  R 

3 2  dBsm2 

Model  1 

NCP  with  b\  —  — 10  dBsm,  q0  =  10° 

Model  2 

NCP  with  b2  =  0  dBsm,  q0  =  10° 

Model  3 

NCP  with  &3  =  0  dBsm,  qo  =  10~3 

Model  4 

NCP  with  b4  =  ±10  dBsm,  q0  =  10° 

The  parameters  for  filters  implemented  in  the  simulation  study  are  shown  in  Table  1 .  The  IMM  filter  was  implemented 
with  four  NCP  models.  Simulations  were  conducted  also  with  an  IMM  with  three  NCP  models  and  one  NCV  model  (the 
NCV  was  to  track  the  slowly  time  varying  RCS).  However,  we  found  slightly  better  results  using  the  four  NCP  models.  In 
the  configuration,  a  central  NCP  model  with  bi  =  0  and  low  process  noise  was  to  provide  the  primary  smoothing  of  the 
mean  (median)  RCS.  Two  NCP  models  with  bi  =  ±10  were  added  to  handle  the  “jumps”  in  the  RCS;  a  larger  family  of 
bias  models  can  be  added,  depending  on  the  magnitude  of  target  RCS  changes  and  the  desired  transition  rate.  A  high  noise 
NCP  model  with  6,;  =  0  was  included  to  handle  conditions  where  the  RCS  was  slowly  varying.  The  following  transition 
matrix  was  used  to  encourage  transitions  to/from  the  central  low-noise  NCP  model. 


0.90 

0.05 

0.05 

0.00 

0.05 

0.75 

0.15 

0.05 

0.01 

0.01 

0.97 

0.01 

0.00 

0.05 

0.05 

0.90 

5.2.  Results 

Simulation  results  for  case  1  are  shown  in  Figure  1.  The  results  are  for  1000  Monte  Carlo  trials.  Plot  (a)  shows  the  mean 
estimated  RCS  over  all  runs  overlayed  on  the  true  mean  RCS.  Plots  (b)  and  (c)  show  close-ups  of  plot  (a)  at  time  50  sec 
and  200  sec.  From  these  plots  we  see  that  the  IMM  filter  has  the  fastest  transition  to  the  steps  in  mean  RCS;  the  median 
filter  appears  to  respond  nearly  as  quickly.  The  a-filter  seems  to  respond  the  slowest,  with  a  “fading  in”  of  the  new  value. 

Figure  Id  shows  the  RMS  value  of  the  RCS  estimate  over  the  1000  trials.  Here  we  see  that  the  IMM  filter  quickly 
hits  an  RCS  floor  value  after  jumps.  This  convergence  is  expected  in  IMM  filters;  inherently,  when  a  small  contribution 
of  the  “non-matching  filter”  is  mixed  into  the  RCS  estimate,  it  “dilutes”  the  accuracy.  Once  can  mitigate  this  by  adjusting 
the  mixing  probabilities  (at  the  expense  of  other  performance  factors),  or  by  choosing  the  “best”  matched  mode  instead  of 
blending  all.  Such  tuning  studies  could  be  pursed  in  the  future.  Meanwhile,  in  Figure  Id  the  median  filter  accuracy  wanes 
after  the  200  sec  and  250  sec  transitions;  during  these  intervals  the  short-window  median  filter  dominates,  and  this  leads 
to  lower  accuracy.  The  a-filter  accuracy  plot  shows  that  it  improves  after  transitions  as  “smoothing”  drives  the  error  down 
during  the  flat  regions. 

Figures  le  through  lg  show  plots  of  the  mode  probabilities  for  the  three  filters  for  case  1.  For  the  IMM,  we  see  that 
the  low-noise  NCP  filter  with  b%  =  0  dominates  most  of  the  time  (red  curve),  which  is  the  desired  response.  At  the  times 
of  the  RCS  jumps,  the  bias  mode  appropriately  peaks  up;  at  time  50  sec,  the  bt  =  10  mode  peaks  up,  at  time  200  sec  the 


b  \  =  —10  mode  peaks  up,  and  at  time  250  sec  the  b  \  =  —10  mode  peaks  up.  We  also  see  that  the  high-noise  NCP  model 
mode  probability  peaks  up  after  these  transitions;  apparently,  this  filter  corrects  for  some  of  the  inaccuracy  in  the  RCS 
estimate  after  the  transition.  Figure  If  shows  the  mode  probabilities  for  the  median  filter.  Here,  the  transitions  are  more 
gradual  compared  to  the  IMM.  Figure  lg  shows  the  mode  probabilities  for  the  a-filter;  they  are  similar  to  the  median  filter, 
but  the  response  is  even  slower.  However,  one  can  see  from  the  plots  that  all  three  filters  properly  “detect”  RCS  jumps  and 
transition  to  the  correct  filter  mode  at  the  transition  point. 

Figure  lh  shows  a  plot  of  the  standard  deviation  of  the  RMS  estimate  generated  by  the  IMM  filter.  Neither  the  median 
filter  nor  the  a-filter  has  the  ability  to  estimate  the  RMS  accuracy  of  the  estimate.  The  plots  shows  that  the  accuracy  estimate 
is  slightly  optimistic  compared  to  the  accuracy  shown  in  Figure  lb.  Nevertheless,  it  correctly  follows  the  error  trend.  The 
availability  of  the  RMS  estimate  is  significant  in  that  this  estimate  could  be  used  in  the  signal  score  (e.g.,  equation  (23)) 
instead  of  a  fixed  parameter.  This  would  enable  finer  tuning  of  the  signal  score,  and  thus  better  differentiation  of  targets 
during  the  data  association  process. 

Plots  for  case  2  are  shown  in  Figure  2.  In  this  example,  the  mean  RCS  varies  with  time.  We  see  that  all  three  filters 
track  the  RCS  reasonably  well.  The  conclusions  about  the  IMM  performance  are  less  clear  here.  In  some  instances  (e.g., 
time  200  sec)  the  IMM  transitions  the  fastest  of  the  three  filters.  In  other  instances  (e.g.,  50  sec),  the  response  is  comparable 
to  the  other  two  filters.  Also,  with  the  time  varying  nature  of  the  mean  RCS,  one  would  expect  the  high-noise  NCP  mode 
(model  2)  to  provide  a  higher  contribution.  A  variety  of  mode  transition  probabilities  and  process  noise  (go)  values  were 
tried,  but  the  low-noise  NCP  filter  with  63  =  0  typically  dominates  in  this  example.  It  is  possible  other  tunings  of  the  IMM 
configuration  would  provide  superior  filter  performance. 

6.  SUMMARY 

In  this  paper,  we  presented  a  methodology  for  using  radar  signal  features  (SNR  and  RCS  measurements  from  a  narrowband 
radar)  for  augmenting  the  track  score  used  in  the  data  association  problem.  Two  main  components  to  the  methodology  were 
the  signal  score  likelihood  ratio  formulation  (Section  3),  and  the  multiple  model  RCS  tracking  technique  to  estimate  the 
mean  RCS  associated  with  an  object  (Section  4).  Signal  scores  were  derived  for  three  different  target  fluctuation  models: 
Swerling  I,  Swerling  III,  and  log-normal.  In  each  of  these  cases,  an  estimate  of  the  mean  target  SNR  (or  RCS)  is  needed. 
The  estimate  is  derived  from  the  RCS  tracking  filter.  Three  different  tracking  filter  options  were  developed:  IMM  filter, 
median  filter,  and  a-filter  (all  using  multiple  models).  A  simulation  study  was  performed  to  compare  the  RCS  tracking 
filters  that  showed  the  following: 

•  All  three  RCS  filters  provide  the  appropriate  tracking  performance:  they  transition  when  large  jumps  occur,  and  they 
provide  smoothing  when  the  RCS  is  relatively  constant. 

•  The  IMM  filter  provides  fast  transition  to  new  values  after  jumps.  The  median  filter  also  typically  provides  a  fast 
transition  after  a  jump.  The  a-filter  typically  transitions  slightly  more  slowly  because  it  “fades  in”  the  values  after  a 
change. 

•  In  terms  of  RMS  accuracy,  the  a-filter  is  superior  because  of  the  “smoothing”  property  during  periods  of  near¬ 
constant  RCS.  The  IMM  filter  can  provide  reasonable  accuracy,  however,  because  of  the  mixing  of  mis-matched 
models  there  is  limit  to  its  accuracy  performance;  techniques  to  improve  the  accuracy  are  available. 

•  The  mode  probabilities  of  all  three  filters  exhibit  the  expected  response  to  jumps  in  the  mean  RCS.  In  some  applica¬ 
tions,  it  could  be  useful  to  monitor  these  mode  probabilities  to  detect  changing  conditions. 

•  The  IMM  filter  is  able  to  generate  an  estimate  of  the  RMS  accuracy  of  the  filter  (from  the  state  covariance).  This  is 
useful  for  computing  the  signal  score  and  may  provide  a  significant  “tuning  capability”  relative  to  the  other  two  filter 
methods. 

In  future  work,  the  RCS  tracking  filter  and  the  signal  score  methodology  will  be  integrated  into  the  tracking  filter  and 
tested  in  closely-spaced  object  scenarios.  This  will  enable  the  value  of  the  feature-aided  tracking  to  be  assessed. 
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Figure  1.  Plots  for  case  1:  (a)  mean  RCS  estimate,  (b)  close  up  of  mean  RCS  at  50  sec,  (c)  close  up  of  mean  RCS  at  200  sec,  (d)  RMS 
of  the  estimate,  (e)  IMM  mode  probabilities,  (f)  median  filter  mode  probabilities,  (g)  a-filter  mode  probabilities,  (h)  IMM  filter  estimate 
of  RMS. 
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Figure  2.  Plots  for  case  2:  (a)  mean  RCS  estimate,  (b)  RMS  of  the  estimate. 
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